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Previous theoretical investigations of colossal magnetoresistance (CMR) materials explain this effect using 
a "clustered" state with preformed ferromagnetic islands that rapidly align their moments with increasing ex- 
ternal magnetic fields. While qualitatively successful, explicit calculations indicate drastically different typical 
resistivity values in two- and three-dimensional lattices, contrary to experimental observations. This conceptual 
bottleneck in the phase- separated CMR scenario is resolved here considering the cooperative nature of the Mn- 
oxide lattice distortions. This induces power-law correlations in the quenched random fields used in toy models 
with phase competition. When these effects are incorporated, resistor-network calculations reveal very similar 
results in two and three dimensions, solving the puzzle. 
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The study of self-organization in transition-metal oxides 
is one of the dominant scientific themes of condensed mat- 
ter physics (CMP). This phenomenon includes the intrinsi- 
cally inhomogeneous states of CMR manganites and un- 
derdoped high temperature superconductors^. In these com- 
pounds, the competition between different ordering tenden- 
cies leads to complexity: their properties change dramatically 
upon the application of relatively small perturbations. In Mn- 
oxides, the cross-fertilization between theoretical and exper- 
imental investigations has been remarkably fruitful, and at 
present the existence of mixed-phase tendencies in the CMR 
regime is widely accepted"^. The emerging CMR picture is 
based on nanoscale clusters of competing phases^-^-^-^. With 
increasing magnetic fields, the clusters with ferromagnetic 
(FM) characteristics rapidly align their moments, leading to 
a percolative insulator-metal transition. Several other com- 
pounds share similar phenomenology, and "clustered states" 
are rapidly emerging as a novel paradigm of CMP«. 

To our knowledge, the only major unresolved issue that 
confronts the phase- separation scenario for the CMR oxides 
concerns the dimensionality dependence of current theoreti- 
cal descriptions. Recent resistor-network calculations in two 
dimensions (2D) reported a colossal MR effect, compati- 
ble with experiments, near the clean-limit first-order FM- 
antiferromagnetic (AF) phase transition. Disorder was fur- 
ther shown to smear the FM-AF transition region into a glassy 
clustered state^. Well-known arguments^ predicts that, in 2D, 
infinitesimal disorder is sufficient to create large coexisting 
clusters of neighboring phases, due to the competition be- 
tween cluster surface effects and the random impurity dis- 
tribution inside a bubble of one phase embedded into an- 
other. However, similar simulations in three dimensions (3D) 
(shown below) do not lead to equally impressive resistivity p 
vs. temperature curves. Within the Imry-Ma reasoning^ the 
critical dimension is 2, and only an unphysically large disor- 
der strength Ac can destabilize the uniform 3D FM phase of 
the Random Field Ising Model (RFIM). Moreover, for A> Ac 
the resulting clusters are not large enough to induce a sub- 
stantial p. Therefore, it is crucial to resolve this incorrect 
dimensional dependence^. Since the phenomenology emerg- 



ing from computer simulations in 2D matches qualitatively 
the experimental results gathered in both 2D and 3D, mixed- 
phase tendencies likely dominate in real materials. More- 
over, recent experiments have unveiled a remarkable insta- 
bility of the CE phase to the introduction of disorder in 3D, 
showing that disorder is by no means irrelevant in real per- 
ovskite manganites This is also compatible with recent 
small-clusters simulations Nevertheless, in the transition 
from realistic models, which cannot be simulated on large 
enough lattices to reach percolation, to the RFIM-like toy 
models that can successfully estimate magnetoresistances&, 
the appearance of an unphysical dimensionality dependence 
suggests that an important conceptual ingredient has been lost. 

In this paper, the dimensionality-dependence puzzle is 
solved. The crucial issue unveiled here is the key relevance 
of cooperative effects for quantitative magnetoresistance stud- 
ies of Mn-oxides. Cooperation introduces correlations in the 
quenched disorder needed to render percolative the clean-limit 
standard FM-AF first-order transition of simple models of 
phase competition. Previous simulations used uncorrelated 
disorder^, and this induced the substantial quantitative dif- 
ferences between 2D and 3D. The disorder discussed here is 
inevitable - and, thus, intrinsic - in the standard chemical- 
doping process widely used to control the hole density, or av- 
erage hopping amplitude, in transition-metal oxides. Replac- 
ing tri- by di-valent ions of different sizes, introduces MnOe 
octahedra distortions that cause local disorder. Once a distor- 
tion is created at a given lattice site, this distortion propagates 
following a power-law decay 1 /r" governed by standard elas- 
ticity mechanisms that suggest a~3 (for recent literature see 
Refs 14 15). This propagation emerges from the cooperative 
nature of the distortions, since adjacent MnOe octahedra share 
an oxygen. This cooperation was already shown to be cru- 
cial for understanding the charge-order states at commensu- 
rate fillings, such as x=0.5. The present effort shows that co- 
operation is also crucial for the understanding of the dimen- 
sionality dependence of CMR simulations. 

The relevance of elastic effects has already been 
emphasized using Ginzburg-Landau free energies. Elas- 
tic compatibility constraints were found to generate texture- 
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FIG. 1 : MC results for the two-orbital DE model with cooperative 
JT phonons, on a 8x8 lattice at low temperature (T). The Hund 
coupling is oo, and the classical spins are assumed ferromagnetically 
aligned for simplicity. Cooperation is included by using the oxygen 
coordinates as d.o.f. Only the Q2 mode is active (not a restrictive 
assumption, since at x=0 and 0.5, Q2 is the most relevant mode^). 
(a) Clean-limit results for A=1.6>Ac=1.4, and x=0.5. The dark and 
open circles indicate positive and negative (Q2), respectively, con- 
comitant with populated orbitals oriented along the x- and ^-axis. 
The (Q2) absolute value, proportional to the dot radius, is related to 
the charge at each site. The previously documented stripes ^'^ "^ are 
clearly observed in the simulation (the small deviations from a per- 
fect arrangement are caused by finite-T effects), (b) Same as (a) 
but with 4 sites (diamonds) where A=0.0. Now the stripe pattern is 
drastically disrupted, showing the high sensitivity of the CO state to 
disorder^ (c) Illustration of the opposite effect as in (a-b): here 
A is subcritical (=1.2) everywhere but in the 4 sites with the largest 
dots where A=2>Ac. A clean-limit study with uniform A=1.2 reveals 
no order, but the inclusion of just 4 sites with A>Ac clearly creates 
short-range stripe order. Here and in (d) the dot area is proportional 
to (Q2). (d) Similar as (c) but for x=0.0. A=0.4 (below Ac=0.5) at 
all sites but the 4 with the largest dots, where A=2.0. The ordered 
plaquette generates charge-ordering on the entire 64-site lattice. 



inducing anisotropic long-range elastic forces, similar to those 
appearing in several transition-metal oxides. Strain effects 
have been discussed in manganite thin films as well^^, and 
they could play a role in stabilizing charge-ordered states 
The effect of long-range Coulomb interactions near first-order 
transitions has also been recently investigated by Yang^^, who 
found an interesting dependence of critical dimensions with 
the interaction range. 

The key relevance of cooperative effects can be dramat- 
ically illustrated via Monte Carlo (MC) simulations of the 
two-orbital double-exchange (DE) model coupled to Jahn- 
Teller (JT) classical phonons. The Hamiltonian and details 
of the simulations have been extensively described in previ- 
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FIG. 2: Top panel: Computer-generated ground- state magnetiza- 
tion {m{T=0)) of the RFIM vs. A, for several values of the disor- 
der correlation exponent. a=2.5 is the predicted 3D critical value 
below which infinitesimal disorder will destroy long-range order. 
Note the dramatic difference among the several as. The inset shows 
{m{T=0)) vs. A for three different lattice sizes and a=3, to illustrate 
finite- size effects. Bottom panel: Snapshots of a typical random field 
distribution hi (top row) and corresponding Ising spin configuration 
(bottom), for three as on a 128^ lattice and A / J=l, to visualize clus- 
ter sizes and shapes. 2D clusters are used to access large linear sizes, 
but results are similar in 3D (see Fig|3j^)i^. 



ous literature and they will not be repeated here. The ex- 
plicit use of oxygen degrees of freedom (d.o.f.) introduces 
cooperation in the distortions. To simplify the calculation and 
allow the study of 64-site clusters, the classical spins were 
frozen in a EM state and the transitions studied here only in- 
volve the charge/orbital d.o.f. This is certainly not restric- 
tive since recent studies have unveiled charge/orbital order- 
disorder transitions at hole-density x=0.5 even with the spins 
in a EM configuration, by varying the electron-phonon cou- 
pling A. The charge-ordered (CO) phase has the same ar- 
rangement of charge and orbitals as the realistic CE state 
Although the MC study in Eig[l]is necessarily restricted to 
2D, this is sufficient to show the key role of oxygen coop- 
eration, illustrating the limitations in previous uncorrelated- 
disorder assumptions. Typical results in the x=0.5 clean-limit, 
and with A larger than the critical value Ac toward a CO- 
state^\ reveal the familiar pattern of charge-diagonal stripes 
with (3x^-r^/3?/^-r^) populated orbitals (Eig[lt). This order 
is dramatically affected when, to simulate quenched disor- 
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der, the value of A is made subcritical in 4 sites of the 64- 
site cluster (just ~6% of the sites, Fig[lJ)). The stripe pat- 
tern disappears and a random-looking distribution of charge 
and orbitals is stabilized, due to the non-local character of 
the disturbance caused by the 4 subcritical sites (compatible 
with recent simulations^ ^•^^). In Fig[l]:, the situation is re- 
versed: the background has a A< Ac and, as a consequence, the 
lowest-energy state is not charge/orbital ordered in the clean 
limit. However, once 4 sites carry A>Ac, a stripe-like pattern 
emerges, affecting most of the lattice. Finally, even at x=0, 
having the 4 sites of a plaquette above Ac - with the rest below 
Ac (see Ref 22)- is sufficient to induce a staggered pattern on 
the entire cluster (Fig [Hi). These realistic-model simulations 
clearly show that cooperation dramatically enhances the role 
of quenched disorder in manganite models. 

The results in Fig [l] indicate that it is inappropriate to use 
the RFIM with uncorrelated disorder to mimic the physics of 
Mn-oxides. If a chemical-doping-induced lattice distortion at 
a Mn-Mn link leads to, e.g., a decrease of the hopping ampli- 
tudes, the neighboring links tend to have a similar reduction 
due to the slow power-law decrease of the elastic distortion. 
As a consequence, a proper RFIM modeling of real mangan- 
ites requires a correlation in the random fields. More formally, 
consider the modified RFIM Hamiltonian 

H = -JY1 S^SJ - A ^ h,Sj/d^ , (1) 

(ij) 

where Si are Ising variables, J is the FM coupling, A is the 
disorder strength, and dij is the distance between lattice sites 
i and j (in practice, dfj was replaced by (1 + df^ )"/^, with 
the same asymptotic behavior but remaining finite at zero dis- 
tance). In this model, a 'random' perturbation hi at site i, in- 
fluences the neighboring dynamical variables sj well beyond 
the usual on-site i=j range, as the analysis in Fig [l] indicates. 
By redefining, hj='^- hi/dfj, the Hamiltonian Eq[l]can be 

cast in the standard form H=—J^^..^ SiSj — A hjSj, but 

now with correlated disorder since (hihj)=l/d'^^~'^ (D = 
lattice dimension). The critical value of a — below which 
the system breaks into domains for infinitesimal A — is 
ac={D/2) + 1, which for D=3 is ac=2.5 (for details, the 
reader should consult previous literature on correlated disor- 
der such as Refl2lL and references therein). The important 
point for our purposes is that correlated disorder can alter the 
critical dimension, and its value can be raised to 3 if a < 2. 5. 
To test these ideas, model EqQhas been studied using algo- 
rithms that allows for the direct calculation of RFIM ground 
states2^ (see Fig|2|i. In agreement with our expectations, there 
is a dramatic qualitative difference between the results ob- 
tained with uncorrelated disorder (a=oo), and those obtained 
using a value of a (a=3) that more realistically mimics the 
elasticity. In particular. Fig |2| shows that the former exhibits 
a large Ac and small clusters, while the latter has large clus- 
ters and a Ac^lO times smaller than the value obtained with 
uncorrelated disorder. 

To further verify these ideas, the formalism presented 
in Ref|^ is here employed. To generate phase competi- 
tion, a simple spin model with competing interactions is 




1 1.2 1.4 1.6 1.8 2 0.6 0.8 1 1.2 1.4 



T/J T/J 




FIG. 3: (a) J' dependence of p (in units of the metallic regions resis- 
tivity po) for the 3D J -J' model with long-range correlated disorder 
(a=3), at W/ J=0.05. (b) p dependence on the disorder strength W, 
at /J=0.23. (c) p in 3D without long-range correlated disorder 
(i.e. at a=oo). Note the small value of p/po compared with (a,b). 
(d) Slice of a 32^ lattice at T/J=1.4, f/J=0.23 and 1^/J=0.05. 
Black and dark-gray represent FM regions with opposite orientations 
of their magnetic moments, and light-gray represents competing state 
(collinear AF) regions. 



used H=-J^^.js^ SiSj-\-j' ^^.^^SiSk, where Si are Ising 
variables, and J (J ) is a nearest-neighbors (next-nearest- 
neighbors) FM (AF) coupling. For small J / J, the T=0 dom- 
inant state is FM, while at large J / J it has collinear AF fea- 
tures (alternating lines of spins up and down). The clean-limit 
critical value is J /J=0.25 (0.5) in 3D (2D). Disorder is in- 
troduced by the replacement J'^J.j^=J -\-Weik at every pla- 
quette diagonal, with Cik being random numbers in [1/2,-1/2] 
spatially correlated as hj . Disorder reduces the values of the 
clean-limit critical temperatures T* to Tc - as extensively dis- 
cussed before^'^ - creating an intermediate T region where 
FM clusters with random moment-orientations are found. A 
grid of resistors can be constructed and the effective cluster 
resistance can be calculated, following standard procedures^. 
Results are in Fig|3l Panel (a) shows the net resistivity p 
vs. T, at several J s, with weak (W<^J) disorder incorpo- 
rated. For all J s, a fairly sharp peak is found between Tc 
and T* for a=3, the exponent that mimics the effect of elastic 
forces. These p-profiles are in good agreement with Mn-oxide 
experiments^^. Figure [S) illustrates the W dependence of the 
results. For sufficiently large W, the clusters are small and p is 
not enhanced at intermediate temperatures. As W is reduced, 
the clusters increase in size and the peak in p develops^^. Fig- 
ure [31: contains p vs. T, parametric with magnetic fields, for 
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FIG. 4: (a) plpo vs. T for the J-J' model (J7J=0.23 and 
VK/J=0.05) with long-range correlated disorder (a=3), in 3D and 
for the magnetic fields indicated. Solid (dashed) lines are results on 
16^ (32^) lattices, (b) Same as (a), but on a 64^ 2D lattice, with 
J7 J=0.68 and 1^/7=0.1. Clearly, now both 2D and 3D results are 
quite similar in magnitude. 

the case of uncorrelated disorder. In agreement with the intro- 
ductory discussion, p here is two orders of magnitude smaller 
than with correlated disorder, illustrating the dramatic differ- 
ences that correlated quenched-disorder causes in the quanti- 
tative results. Finally, dominant MC configuration snapshots 
(FiglSji) reveal an intricate cluster arrangement in the T re- 
gion of interest if a=3 — far from the uniformly polarized 
state at the same W if a=OQ — intuitively justifying the ob- 
served high p values2i&. The proximity of a=3 to ac=2.5, and 
the anticipated further effective reduction of a if Coulombic 
disorder effects were incorporated, lead us to believe that the 
sub-/im clusters reported by Uehara et ah (RefU) could in- 



deed be intrinsic to Mn-oxides. 

Although the previous results clearly show that disorder- 
correlation effects are crucial, only a magnetoresistance esti- 
mation can clarify its role in the CMR effect. For this pur- 
pose. Fig. |4| contains a 3D/2D resistor-network calculation of 
p, with correlated disorder, varying the magnetic field. The 
peak in the resistivity - in the region between Tq and T* - 
is quite clear and has a similar value (just a factor 3 of differ- 
ence) in both 3D and 2D. This is a considerable improvement 
over results with uncorrelated disorder, where the peak resis- 
tivity ratios between 2D and 3D are as high as 200 or more 
(see Ref 6 and Fig|3l:). The effect of small magnetic fields 
- that rotate large preformed FM clusters - is now strong in 
both dimensions of interest, and colossal MR ratios are ob- 
tained with minimal tuning of couplings, as shown in Fig|4| 

Summarizing, here the key role of cooperative effects in 
the theoretical description of the CMR effect has been un- 
veiled. Cooperation induces correlation in the disorder needed 
to transform a first-order FM-AF transition into a percolative 
process. Explicit calculations in toy models for the CMR phe- 
nomenon show that the critical dimension of the system is al- 
tered by disorder correlation, and when elasticity effects are 
included magnetoresistance ratios of comparable magnitude 
are obtained in 2D and 3D. These results remove a concep- 
tual roadblock of previous phase- separation-based theoretical 
studies of manganites by demonstrating the importance of cor- 
related disorder induced by cooperative strain effects, reaf- 
firming the relevance of clustered states in the description of 
transition-metal oxides. 
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The p reduction as the critical J =0.25 is reached in (a) or as 
W is reduced in (b) is a size effect, which occurs when intrinsic 
coexisting-cluster sizes become comparable with MC lattices. 



